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Abstract 

We perform a combined next to leading order analysis to both inclusive and 
semi-inclusive polarized deep inelastic scattering data, assessing the impact of new 
data in the extraction of the sea quark and gluon spin dependent densities. In 
particular, we find that semi-inclusive data show a clear preference for a positive 
Au distribution, although are almost insensitive to different alternatives to Ad. 
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1 Introduction 



Over the last decade, there has been a remarkable improvement in our 
knowledge of the spin structure of the nucleon at partonic level, driven both 
by the increase in precision and volume of data on spin dependent structure 
functions, and also on theoretical and phenomenological work around the ever 
increasing wealth of experimental evidence |l[ . In the process, the extraction 
of spin dependent parton densities has evolved from mere extensions of simple 
scale independent SU(6)-inspired models, to consistent next to leading order 
(NLO) QCD global analysis based on a much more reduced set of model 
dependent hypotheses. 

In spite of this significant progress, until very recently it was almost un- 
avoidable to make rather crude assumptions about the way in which sea 
quarks are polarized, or more precisely, how polarization is distributed among 
the different sea quark flavours. The main reasons for this being the inability 
of totally inclusive data to discriminate between quarks and antiquarks, and 
the lack of enough statistics in the pioneering polarized semi-inclusive ex- 
periments. Because of these circumstances, up to now even the most refined 
global analyses have simply assumed polarization to be evenly distributed 
among the sea flavours, or in the most ambitious attempts, they have tried 
to discriminate between light and heavier flavours, but always in compliance 
with SU(2) isospin symmetry assumptions @, |3|, |], || . 

Motivated by the growing evidence supporting an SU(2) symmetry break- 
ing between unpolarized sea quark distributions, more recently there has 
been a great deal of activity concerning sea quarks polarization, including 
theoretical and phenomenological models advocating different degrees and 
realizations of SU(2) symmetry breakdowns |J, and even attempts to find 
evidence for them in the data f7|. 

Semi-inclusive asymmetries are sensitive to contributions coming from 
specific combinations of quarks of different flavour and nature, so in principle 
they can shed some light on the sea quark polarization issue. The theoretical 
framework needed to deal with them consistently up to NLO accuracy is 
well known ||, and indeed this kind of data was successfully included in 
a NLO global analysis some time ago |§. However, the statistical weight 
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of the data available at that time H hindered a significant impact on the 



parton distribution extraction. Since then, the Hermes Collaboration |T0 
have presented new semi-inclusive data with greater precision, circumstance 
that justifies a new analysis taking into account also the new additions to 
the inclusive data set. 

In this paper we perform a combined next to leading order analysis of 
both inclusive and semi-inclusive polarized deep inelastic scattering data, 
releasing the usual assumptions about flavour symmetry relations between 
the polarization of both valence and sea quarks, and we assess the role of 
semi-inclusive data in the experimental evaluation of these relations. 

Performing this analysis we have found that semi-inclusive data consti- 
tute a clear constraint on the polarization of 'up' antiquarks, information 
which combined with that coming from inclusive measurements, fix the po- 
larization of valence quarks of the same flavour. However, 'down' antiquarks 
polarization remains loosely bounded allowing a large uncertainty on 'down' 
valence quark densities. As in previous analyses, the data imply a weak 
constraint on the polarized gluon density. 

This paper is organized as follows: In section 2 we review the framework 
for global analyses of inclusive and semi-inclusive data, and present our con- 
ventions. In sections 3 and 4 we present the result for the fits to inclusive 
and combined (inclusive plus semi-inclusive) data, respectively. Finally, in 
section 5, we give our conclusions and in the appendix we present the pa- 
rameters for the fits. 



2 Conventions and Framework 



Throughout the present analysis we adopt the MS factorization scheme, 
as implemented in our previous analysis 0, for writing inclusive and semi- 
inclusive asymmetries in terms of parton distributions up to NLO accuracy. 
For the inclusive asymmetries we take as usual 

A N (r O 2 ) ~ 9i{x,Q 2 ) m 



2 



where the spin dependent structure function is given by 

ni{*, Q 2 ) = \y, 4 { q 2 ) + ® ^ + a q, ® A ^]} > ( 2 ) 

and in the case of three flavours can be alternatively written as 

^,Q 2 ) = (±^ s + l 6 ^ s + 1 as) ® (l + gAC ? ) 

+ £e^A<7®AC fl . (3) 

g 

In equation (3) the ± sign corresponds to either proton (N = p) or neutron 
(N = n) targets and 

(Au + Am) - (Ad + Ad) 
(Am + Am) + (Ad + Ad) - 2 (As + As) 
(Am + Aw) + (Ad + Ad) + (As + As) . (4) 

The first two (non-singlet) distributions evolve in Q 2 independently, whereas 
the gluon and the singlet distributions are coupled in the evolution. 

It is worth noticing that whereas Aq^ s can be obtained directly from 
data on the difference between the proton and the neutron spin dependent 
structure function g\ — g™ , the three remaining distributions that appear 
in equation (|3]) have to be obtained somewhat indirectly, i.e. exploiting 
their different scale dependence. Unfortunately, the limited amount of data 
available, and the rather restricted range in x and Q 2 covered, do not allow to 
do a precise determination of all the distributions. Nevertheless it should be 
kept in mind that even with an unbounded set of inclusive data of unlimited 
precision, it would only be possible to determine combinations of Ag + Aq = 
Aq v + 2Aq (with q = u,d,s), but never Aq v and Aq separately, as can be 
seen from equation (|j); therefore, the data is insensitive to the symmetry 
breaking in the sea sector.^] 

3 In this sense, claims stating that 'the data tolerate strong flavor symmetry violations 
in the sea sector' (as in reference ||) are just rather void statements, since of course it is 
always possible to modify Aq in such a way that Ag + Aq remains unchanged. 
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To be able to determine both Aq v and Aq one has to introduce some 
other independent input in the analysis, such as weak structure functions 
and Drell-Yan production, for which the theoretical framework is well known 



11, 12] but unfortunately there is no data yet, or semi-inclusive deep inelastic 



scattering (SIDIS) asymmetries. 

In order to perform the extraction of polarized valence and sea quark 
densities, we include SIDIS asymmetries measured by Hermes and SMC, 
which are given by 

A ' ^Q) z - SzdzF r(x,z,Q*y (5) 

There, the semi-inclusive analog of gf is given by 

g? h (x, z, Q 2 ) = £ A {Aftfo Q 2 )D h Jz, Q 2 ) 



1,1 

^ aJQ 2 ),. . _ _ fc . 



^4^-[A qi ® AQ, ® D h q , + Aq t ® AC ig ® D h g + Ag <g> AC gj <g> D%\ 



+ AM^x, z, Q 2 ) + ^?S[AM* ® AC; + A< ® AC g ]| , (6) 

and Z denotes the kinematical region covered by final state hadrons [z^ = 
P ■ h/P ■ q > 0.2 in |9|, [UJ). The first line in equation (^) represents the 
familiar LO contribution to semi-inclusive processes, the second accounts for 
current fragmentation up to NLO, and the third is associated with target 
fragmentation and can be neglected for the kinematical region explored by 



Hermes and SMC The full expressions for the coefficients in equations 



(2) and (6) in the MS can be found in || and references therein. 

SIDIS data are available for different kind of targets and species of hadrons 
identified in the final state (p,D,He and ir^, respectively). Since frag- 
mentation functions depend strongly on the flavour of the originating par- 
ton, and on the charge of the hadron, it is possible to construct alternative 
combinations of parton distributions which are independent from the ones 
contributing to the inclusive asymmetry. In this way, in principle it possible 
to do a full flavour decomposition for the polarized parton distributions. 

Fragmentation functions are obtained from fits to e + e~ — > h data and, 
similarly to what happens in inclusive DIS, those measurements allow only 
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the determination of the the combination + D^. In order to construct 
the full set of fragmentation functions needed to analyze the SIDIS data for 
h + and h~, one has to do some extra assumptions. In the present analysis 
we take into account only contributions from 7r's and i^'s, and for both 
mesons we define a Q 2 and flavour independent function which allows 
to separate between the 'favored' and 'unfavored' fragmentation functions 
for each flavour and final state hadron as 



and 
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(8) 



In agreement with measurements of unpolarized SIDIS data |fI5 
eterize r] D as 



, we param- 



Vd(z) 



l + z 



(9) 



The usual Field & Feynman parameterization for r/o [|14| corresponds to 
equation (||) with a set to 1. In our case we will vary the parameter a around 
the default value of 1, as a way to quantify the uncertainties introduced by 
the assumptions on the fragmentation functions in the extracted polarized 
parton distributions. 

Regarding the normalization of equations ([!]) and (|^), is customary to 
write the unpolarized structure functions Ff(x,Q 2 ) in terms of the more 

Therefore the Q 2 dependence of the 'favored' and 'unfavored' fragmentation functions 
is given, in our approximation, by the one of the total distributions fitted from e + e~ data 
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familiar F 2 N (x, Q 



p ( x Q 2 ) f 2 (x } Q 2 )[i + 1 2 (x } Q 2 )} (1Q) 

2x [l + ^g 2 )]' 1 j 



where 7 2 (x,Q 2 ) = 4M 2 x 2 /Q 2 , R(x,Q 2 ) is the ratio between the longitudi- 
nal and transverse photoabsortion cross sections, and we have dropped the 
nucleon target label N. Therefore, the asymmetry for inclusive data is given 
by 

a n*\- 9i(x,Q 2 )2x[l + R(x,Q 2 )] 
Al[X ^ } F 2 (x,Q 2 ) [1 + 7 2 (*,Q 2 )] ' 
and similarly for the semi- inclusive one. 

The information on the polarized parton distributions, whose extraction 
is the goal of the global fit, is therefore fully contained in the polarized 
structure function g±(x, Q 2 ), which can be reconstructed using the known 
(measured or computed) values of Ax(x, Q 2 ), F 2 (x,Q 2 ) and R(x,Q 2 ). 

Of course, the Ai(x,Q 2 ) values are obtained from measurements, but 
there is still some freedom in taking F 2 and R, i.e. either from experimental 
data or from pQCD calculations. Since both quantities are strongly affected 
by higher twist contributions and unknown higher order corrections in the 
measured range of x and Q 2 , but only a leading twist LO/NLO expression 
for gx is used in equation (O) under the assumption that higher twist effects 
may cancel when taking the ratio, in the past it has been considered more 
consistent from the perturbation theory point of view to use the pQCD result 
for both F 2 and R. 

On the other hand, it is worth noticing that the measured quantity in 
polarized DIS is Au, related to A% by 

A\\ = D(A 1 + r]A 2 ) ~ DAx, (12) 

where the depolarization factor D is obtained from the measured value of 
R. Furthermore, it turns out that the effect of using either R exp or Rqcd 
can represent a difference up to 20% (40%) at NLO (LO) for the smallest 
values of measured x, and still persists, although becomes smaller, even for 
values of x and Q 2 where higher twist effects are expected to be small (like 



6 



for SMC data). An analogous effect is expected in F2, however in this case 
the differences are found to be smaller. 



In our analysis we have decided to use the experimental value of R from 
reference in order to be consistent with the extraction of A\ from A\\, 
but still construct F 2 using the pQCD expressions that can be obtained from 
the unpolarized analog of equations (g) and © and the GRV98 JlTj set of 



unpolarized parton distributions. The choice for a perturbative F 2 , is forced 
by the fact that our fit is also performed to semi-inclusive data, and thus we 
need consistency between the normalization of inclusive and semi-inclusive 
asymmetries. For the time being no experimental values of F£ are available 
for the normalization of the latter, at least for the same bins of x and Q 2 
measured by SMC and Hermes. Regarding R h exp , we have just assumed that 
the same value applies for both inclusive and semi-inclusive data. 

Clearly, the extraction of R h and F£ from unpolarized semi-inclusive data 
will be a great advantage in the future in order to perform a fully consistent 
analysis. Nevertheless, as we have stated before, there are still some other 
uncertainties in the semi-inclusive case which are expected to hinder the ones 
coming from eventual differences between R and R h . 

We also implement a different strategy to that used in || for parameter- 
izing and extracting parton densities. The procedure consists in two stages. 
The first one, in which totally inclusive data is used to fix the net polar- 
ization carried by each quark flavour, but without discriminating between 
quarks and antiquarks, as this distinction is beyond the reach of this kind 
of data. For u and d quarks plus antiquarks densities at the initial scale 
Qq = 0.5 GeV 2 we propose 

x aq (l- xY q (l + ^ q x 5q ) 
X{Aq + Al!) = N<1 B(a q + 1,^ + 1) + 7,5(^ + ^ + 1,^ + 1)' (13) 

while for strange quarks plus antiquarks we use 

x(As + As) = 2N \ } (14) 
B{a s + l,p s + l) 

with a similar parametric form for gluons 

xA g = N q f 1 J -. (15) 

y g B(a g + l,{3 g + l) 1 ; 
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The first moments of the quark densities 5q (N q ) are often related to 
the hyperon beta decay constants F and D through the SU(3) symmetry 
relations 

5u + Su - 8d - 5d = N u - N d = F + D = 1.2573 (16) 

5u + 5u + Sd + 5d - 2(5s + 5s) = N u + N d - AN S = 3F - D = 0.579. (17) 

Under such an assumption, the previous equations would strongly constrain 
the normalization of the quark densities. However, as we are interested also 
in testing flavour symmetry at different levels, we leave aside that strong as- 
sumption and relax the symmetry relations introducing two parameters, e#j 
and e S u(3) respectively. These parameters account quantitatively for even- 
tual departures from flavour symmetry considerations, including also some 
uncertainties on the low-x behavior, and higher order corrections, 

N u -N d = (F + D)(l + e Bj ) (18) 

N u + N d - AN S = (3F -D)(l + e Sf/(3) ), (19) 

and we take them as a measure of the degree of fulfillment of the Bjorken 
sum rule [|19| and the SU(3) symmetry. 

Equations ( jI5[) and (|T9| ) allow to write the normalization of the three 
quark flavours in terms of N s , 6b j, and esu(3)- Notice that no constraints have 
been imposed on the breaking parameters since we expect them to be fixed by 
data. The remaining parameters are constrained in such a way that positivity 
with respect to GRV98 parton distributions is fulfilled. This is particularly 
relevant at large x, and since no polarized data is available in that kinematical 
region, we directly fix the parameters (3 U = 3.2 (3.05), (3 d = 4.05 (3.77) and 
P g = 6 (6) for the NLO (LO) sets in agreement with GRV98. Consistently 
with the choice for the unpolarized parton distributions, we use the values 
of Aqcd given in Ref.|T^ to compute a s at LO and NLO. 

In this way, in the first stage we obtain from the inclusive data set all the 
information available there without any unnecessary symmetry assumptions. 
The second stage of the procedure incorporates semi-inclusive data in order to 
extract antiquark densities, and consequently to fix valence quark densities, 
relaying on the possibility to discriminate between quarks and antiquarks 
given by positive and negatively charged hadron production asymmetries. 
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As antiquark densities we take 



*** =J W l,ft + lr (20) 

for u and d quarks, and we assume s = s since the possibility of discrimination 
in the s sector is beyond the precision of the data and j3g = (3 S . 

The data sets analyzed in both stages include only points with Q 2 > 1 
GeV 2 listed in Table 1, and totaling 137, 118, and 34 points, from proton, 
deuteron, and helium targets respectively, in the inclusive stage, plus 42, 24, 
and 18, from proton, deuteron, and helium targets respectively, in the sec- 
ond stage. In order to take into account correlations between inclusive and 
semi-inclusive measurements, we only analyze the inclusive data for SMC and 
Hermes for 'averaged' bins. In the semi- inclusive case we take into account 
in the fit only the most precise data concerning the production of charged ± 
hadrons (without identifying pions, kaons, or other particles individually). 



Collaboration 


Target 


Final state 


# points 


Rcfs. 


EMC 


proton 


inclusive 


10 




21 




SMC 


proton, deuteron 


inclusive 


12, 12 




22 




E-143 


proton, deuteron 


inclusive 


82, 82 




23 




E-155 


proton, deuteron 


inclusive 


24, 24 




24 




Hermes 


proton, helium 


inclusive 


9, 9 




10 




E-142 


helium 


inclusive 


8 




25 




E-154 


helium 


inclusive 


17 




23 




SMC 


proton, deuteron 


h+, hr 


24, 24 




Hermes 


proton, helium 


h+, h- 


18, 18 




1Q| 



Table 1: Inclusive and Semi-inclusive Data used in the fit. 



3 Inclusive Data 



In this section we present the results coming from the analysis of the 
inclusive data set. Taking into account that in previous analyses, even those 
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which have used most recent inclusive data, the gluon densities have been 
found to be rather poorly constrained, we extract three different sets with 
the gluon first moment constrained to be lower than 0.4 and labeled as 'i', 
bounded between 0.4 and 0.7 and denoted as 'ii', and greater that 0.7 labeled 
as 'iii', respectively, at the initial scale. Doing this, we aim to explore different 
regions for gluon polarization and, at the same time, to provide different 
scenarios for predictions in future experiments. As we will show, data prefer 
gluon distributions with first moments around those values but does not 
really constrain them in a precise way. 




The resulting densities from the LO and NLO analysis, and the corre- 
sponding relevant features at 10 GeV 2 , are shown in Figure 1 and Table 2, 
respectively. The corresponding parameters are given as an appendix (Table 
A3). 

Even though the precision attained by this kind of fits is remarkable good, 
as shown in Figure 2, there seems to be a slight systematic discrepancy be- 
tween Hermes data and the results of the remaining experiments. For these 
reason, in the fitting procedure we have applied to Hermes data a normal- 
ization factor which we allow to depart from unity up to a 5% (it saturates 
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at the upper limit), without almost any consequences in the obtained parton 
distributions, but with a significative improvement in the y 2 values. 



10' 



10 



10 

X 

rT 10 

a 

O. — 
(50 j 



10 



10 




x = 0.005 (f=6400) 
— x = 0.008 (f=3200) 

x = 0.015 (f=i600) 
x = 0.025 (f=800) 
x = 0.035 (f=400) 

x = 0.05 (f=i80) 



1 EMC 
o E155 
' E143 
• SMC 

A HERMES 




x = 0.08 (f=80) 
-*- x = 0.125 ( f=40) 

= 0.150 (f=20) 

0.185 (6=10) 

^ »- x = 0.25(f=5) 

- x = 0.30 (1=2.5) : 

+ *- x = 0.35 (i=i.3) " 

x = 0.416 (f=o.8j 
4— x = 0.468 (1=0.5) 

-* X = 0.5(1=0.25) 

x = 0.75 (1=1.0) : 
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Figure 2: gi data against NLO inclusive set i) 



As can be seen in Figure 1 the three sets have rather large differences in 
the gluon distribution (with a very mild preference for set T at NLO and 
'ii' at LO), with close results for net polarization of each quark flavour. The 
quality of the fits are comparable, in fact, it is hard to distinguish between 
the inclusive asymmetries estimates coming from these sets. As it is shown 
in Table 2, the differences in x 2 are rather small, and allowing variations 
within just two units, it is possible to go from an almost negligible polarized 
gluon distribution to a highly polarized one. The reason for this can be easily 
understood from Figure 2, where all proton data for g\ is shown as a function 



11 



of Q 2 . DIS data, in principle, constrain the gluon distribution through the 
evolution of g±, but unfortunately, available data do not show a real lever-arm 
in Q 2 , particularly at lower values of x where the gluon density is expected 
to drive the evolution. 

Even though available data indicates that Sg(Q 2 = 10) ~ 1 ± 1 (see 
Altarelli et. al. in reference [|| for a more exhaustive analysis), we clearly 



have to wait for less inclusive data as in jet [pTfl , prompt photon p§| , and 
heavy quark |29| production to be able to determine the distribution with a 
reasonable precision. 

In Figure 1 we have included also u, d, and s distributions obtained in 
our previous analysis (DSS fit 3)||. As can be observed, the differences in 
the u and d distributions are rather small, and although there have some 
additions of new data in the analysis, they are mainly originated in the use 
of the experimental value of R instead of the perturbative one as it was done 
in our previous analysis. This is particularly manifest in the LO case, where 
the perturbative approach yields R = 0. 

It is also worth remarking that the fits prefer d distributions steeper at 
small x than the u ones, i.e. < a u . Since there is almost no neutron data 
able to constrain Ad + Ad at x < 0.01, this effect should rather be considered 
as a result of the extrapolation of the distributions than as an experimental 
evidence of g™ being steeper than g\. Clearly, more precise neutron and 
deuteron data at small x would help to clarify this feature. 

Finally, the most prominent difference between the results of the present 
analysis and those from DSS, is the one found between the strange spin 
dependent distributions. Even though their first moments are both negative 
and rather close in magnitude, the shape of the distributions differ, mainly 
because As is not bounded to obey the SU(3) symmetric relations within the 
sea, constraint that was applied in the DSS analysis [[|. This difference give 
also an idea of the large uncertainties on this distribution, which is clearly 
not well constrained by the available data. 

Regarding the symmetry breaking for inclusive data, it is worth noticing 
that for the three NLO sets the Bjorken Sum Rule breaking parameter €bj is 
found to be of the order of 1 %, while the SU(3) breaking parameter esu(3) 
varies around a few percent, and can not be determined with good accuracy 
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from the available data. LO sets show similar features although with slightly 
higher x 2 , and larger symmetry breaking parameters. The larger value of esj 
at LO is mostly related to order a s corrections included at NLO but absent 
there. 



set 


x] 


e Bj 


e SU(3) 


5u + 5u 


Sd + Sd 


5s 


Sg 


5E 


NLO ii 
iii 


228.66 


-0.015 


-0.003 


0.780 


-0.458 


-0.064 


0.75 


0.194 


228.69 


-0.018 


0.037 


0.776 


-0.458 


-0.071 


1.07 


0.175 


230.33 


-0.002 


-0.063 


0.767 


-0.486 


-0.066 


1.75 


0.15 


LO ii 
iii 


237.55 


-0.126 


0.016 


0.713 


-0.385 


-0.065 


0.88 


0.199 


236.10 


-0.126 


0.018 


0.708 


-0.39 


-0.068 


1.4 


0.181 


236.44 


-0.126 


0.009 


0.702 


-0.396 


-0.069 


1.68 


0.169 



Table 2: First moments for LO and NLO distributions at Q 2 = 10 GeV 2 

In Table 2, we also show at Q 2 = 10 GeV 2 the results for the first moment 
of each flavour distribution and the singlet <5£, which is found to be between 
0.2 and 0.15 and correlated to the value of 5G. 



set 




r? 


™ 
1 i 


NLO ii 
iii 


0.191 


0.130 


-0.061 


0.191 


0.129 


-0.062 


0.194 


0.126 


-0.067 


LO ii 
iii 


0.183 


0.130 


-0.053 


0.183 


0.128 


-0.055 


0.183 


0.126 


-0.057 



Table 3: First moment of g\~ n ^ n at Q 2 = 10 GeV 2 

Finally, in Table 3 we present the LO and NLO results for the first mo- 
ments Fi' n = J dxg[ ,n and the Bjorken sum rule evaluated at Q 2 = 10 GeV 2 . 
The differences in the results obtained with the three sets are mostly due to 
the small x extrapolation, in the case of the proton, and the lack of enough 
precise data, for the neutron targets. 
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4 Combined data 



Having obtained all the information available from the inclusive data, we 
proceed to analyze also the antiquark densities, exploiting the semi-inclusive 
data. The constraining power of these asymmetries can be schemed in the 
following way. The production of positively (negatively) charged hadrons 
from proton targets is dominated by the convolution of Au {Ad) densities 
and 'favored' fragmentation functions Df, weighted with a factor of four 
(one) due to the electric charge, 

2<?r p +H ~ ^Au®Dl {2) + l -Ad®Dl (l) 

+^Au <g> Dl {l) + ^Ad ® Dl {2) + ^ As <g> D 3 , (21) 

where we can neglect the suppressed strange quark contributions, and in 
order to simplify the discussion, we have taken into account only final state 
7r's and within a LO expression. 

Factoring out the q + q distributions which are almost completely fixed 
by the inclusive data we can write 

2^ +( - ) ~ ±(A U + Au)®D* {2) + ±(Ad + Act)®D% 1) 

+i(Ad - 4AU) ® (Dl {2) - (22) 

where it is clear that the unconstrained component of both asymmetries is 
proportional to the same factor (Ad — 4Au) but with opposite signs, and, 
regardless the values of D\ and D2, proton asymmetries are significantly more 
sensitive to Am due to the difference in the electric charge factors. 

Semi-inclusive asymmetries from neutron targets are proportional to a 
complementary combination of Ad and Au, and therefore are strongly depen- 
dent on Ad. However, in practice, these last asymmetries have little impact 
in the analysis since they are heavily suppressed by the normalization coming 
from deuterium and helium unpolarized structure functions. Consequently, 
although Au can be more easily pinned down, Ad densities are not very well 
constrained by present data. Indeed we observe this feature in the global 
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fits: the impact of Ad in terms of x 2 is very small and a large modification 
in its distribution can be easily balanced by a small modification of Au. 

Taking into account the previous discussion, we parameterize Au as in 
equation (14) (fixing a„ = a-g) but we restrict Ad to two somewhat extreme 
alternatives: one, labeled as '+', in which = N-j, and another, labeled 
as '— ', in which = —Nj. These scenarios are motivated by expectations 
from different models that have been discussed in the literature. The first 
one (Ad = Au at Q§) corresponds to the usual assumption about SU(2) 
symmetry between light sea quarks. The second alternative (Ad = —Au at 
Qq) breaks the SU(2) symmetry at sea level along the lines of the statistical 
and chiral quark soliton models |EJ. 

The results of allowing both the '+' and the ' — ' alternatives in the ex- 
tended analysis that incorporates SIDIS data are presented in Table 4, where 
we show only the results associated to Sq and Sq v since those for the remain- 
ing parton densities are almost equal to the 'inclusive' ones of Table 2. This 
similarity can be easily understood taking into account the greater precision 
and the amount of these data. The parameters for the NLO and LO fits are 
presented in Tables Al and A2 in the appendix, respectively. 



set 


X 2 T 


Suy 


Sdy 


5u 


6d 


NLO 


i+ 


330.52 


0.649 


-0.576 


0.063 


0.061 




330.00 


0.686 


-0.357 


0.045 


-0.049 


ii+ 


330.98 


0.656 


-0.555 


0.058 


0.056 


ii— 


330.22 


0.684 


-0.357 


0.043 


-0.046 


iii+ 


332.31 


0.652 


-0.583 


0.055 


0.053 


m— 


331.64 


0.679 


-0.392 


0.041 


-0.044 


LO 


i+ 


336.24 


0.599 


-0.493 


0.060 


0.060 




335.54 


0.629 


-0.291 


0.043 


-0.043 


ii+ 


334.19 


0.592 


-0.498 


0.058 


0.058 


ii— 


334.03 


0.624 


-0.295 


0.042 


-0.042 


iii+ 
iii— 


334.26 


0.591 


-0.490 


0.056 


0.056 


333.69 


0.624 


-0.298 


0.041 


-0.041 



Table 4: \ 2 of the LO and NLO fits and the first moment of the distributions obtained 
from the combined analysis at Q 2 = 10 GeV 2 . 

In Figure 3 we show the corresponding parton densities (only for sets 
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i+ and i— , as the other ones are very similar), and in Figures 4 and 5, the 
comparison with semi-inclusive and inclusive data, respectively. The lines in 
Figures 4 and 5 interpolate, as usual, the results of the fit at the measured 
value of Q 2 of each data point, and for E-143 we show only the averaged data. 
For the experiments that provide only the values of g\, we have computed 
the corresponding asymmetries with the procedure discussed above. 
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- 1 1 r 
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V 
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Figure 3: NLO polarized densities from fits to combined data at Q 2 — 10 GeV 2 

In Figure 4 we include also an estimate for the semi-inclusive asymmetries 
computed using the parton distributions obtained in the analysis of inclusive 
data and assuming SU(3) symmetry of the sea, i.e. Am = Ad = As. As it 
can be noticed, the differences are significative in the case of Hermes data 
from proton targets, which completely dominates the fit of semi-inclusive 
data, whereas no distinction can be observed in the almost vanishing He 
asymmetries. 

Furthermore, as can be expected from equation fl22[) the effect of a positive 
Au distribution is to decrease (increase) the asymmetry for the production 
of positively (negatively) charged hadrons from proton targets with respect 
to the case of the SU(3) symmetric result (which of course corresponds to a 
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negative Au), as required by the experimental data. 




10 " 2 10 _1 10 _1 1 

X X 

Figure 4: Semi-Inclusive asymmetries against NLO sets 

As we mentioned, in the computation of the total x 2 values we properly 
take into account the correlations between SIDIS and inclusive data, finding 
that its main effect is just to increase slightly the total x 2 of the fits. 

Notice that the six sets of parton densities in Table 4 have positive Au 
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-2 -1 -2 -1 

10 10 1 10 10 

X 

Figure 5: Inclusive asymmetries against NLO sets 

densities, however their sizes are still correlated to the choice for Ad. This is 
in agreement with the discussion below equation ( p^) , since the correlation 
between Am and Ad is such that (Ad — AAu) remains roughly the same for 
all sets. 

It should be noticed also that even in the case of S'?7(2)-symmetric alter- 
natives (labeled as '+'), there are small differences between the first moment 
of the Am and Ad at Q 2 = 10 GeV 2 , notwithstanding they are set to be equal 
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at the initial scale Qq. The difference comes from the fact that these dis- 
tributions evolve differently at NLO due to both the appearance of different 
non-singlet quark splitting functions and the non-zero value of Au v — Ad v . 
Therefore the SU(2) symmetry in the light sea sector is perturbatively bro- 
ken, as it happens in the unpolarized case, but the size of this breakdown 
is small compared to what is required by the SIDIS data. Notice that this 
effect is absent in the LO case, where there is only one qq splitting function. 

It is worthwhile mentioning that the slight differences in the \ 2 between 
'+' and '— ' sets are comparable to the uncertainties of the fitting procedure, 
such as the use of slightly different parameterizations for the initial distri- 
butions, the arbitrariness in choosing the initial scale Qq, or other external 
constraints. For this reason, we can conclude that Ad densities are not well 
constrained by present data, and the slight differences in \ 2 fail to justify a 
definite statement about isospin symmetry at sea quark level. At variance 
with SU{2), SU(3) symmetry in the sea seems to be clearly broken in the 
sense that fits prefer Am and As with opposite polarization. 

In order to try to quantify uncertainties in the spin dependent distri- 
butions, specifically, those introduced by the assumptions done about the 
fragmentation functions, we have varied the parameter a of equation (|9|), 
which was fixed to 1 in the previous fits, from 0.75 to 1.25. Doing this, we 
find that this variation produces only negligible changes in the Au distri- 
bution, which is explained by a large cancellation of its effect in the rate 
of polarized and unpolarized SIDIS structure functions. Of course, this by 
no means constitute an exhaustive study of the related uncertainties, but it 
makes us more confident about the sensitivity of our results to some unavoid- 
able assumption. 

Although for the time being we do not intend to perform a comprehensive 
analysis of other uncertainties involved in the the extraction of the parton 
distributions presented here, problem which is far from being solved even 
in the unpolarized case we find reasonable to claim that present data 
give a clear signal of a positive Am distribution. Actually, we have found to 
be impossible to produce parton densities constraining All to be negative, 
unless we allow an increase in the \ 2 value of the order of 5% (about 10-20 
units), showing the sensitivity of the semi-inclusive data. 
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Noticeably, most of the available theoretical calculations M predict a 
positive Am distribution. In Figure 6 we show our LO results for Am from 
fits 'i+' and 'i— ' compared to the prediction of Gliick and Reya, and the 
one of Bhalearo || at Q 2 = 1 GeV 2 . As can be seen, there is a reasonable 
agreement between those predictions and our results in the kinematical region 
covered by Hermes semi-inclusive data (i.e. 0.4 > x > 0.03), at least within 
the expected uncertainties on the extracted distributions. Of course, there 
is no experimental evidence yet at smaller values of x so our results just 
provide an extrapolation in that region. The same models predict also a 
negative Ad, however as stated before, global fits still do not discriminate a 
preferred alternative to compare with. 



0.03 

0.025 

x Au 

0.02 

0.015 [ 

0.01 
0.005 



set i (+/-) (LO) 
Gliick and Reya 
Bhalerao 




Figure 6: xAu (set i) at Q 2 = 1 GeV 2 compared to the predictions in Ref. || 

At this point it is worth comparing our results to the ones reported by 
the Hermes Collaboration in reference [nj. There, no definitive sign for Au 
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is found. Besides the fact that we take into account the correct Q 2 evolution 
of the distributions both up to LO and NLO accuracy, whereas Hermes uses 
the naive parton model expressions assuming Q 2 independence of the asym- 
metries, the main difference comes from the fact that we perform the fit to 
all available inclusive and semi-inclusive data (not only the one coming from 
Hermes), and that we do not impose any constraint on Aq distributions. In 
the Hermes analysis sea quark distributions are constrained to have the same 
sign, and taking into account that the strange quark distribution As is found 
to be negative already from the inclusive data, this constraint introduce a 
strong bias against the £77(3) breaking in the found in our analysis. 

5 Conclusions 

The main conclusion of this analysis is that semi-inclusive spin dependent 
data available today, are not only in perfect agreement with inclusive data 
with regard to the extraction of polarized parton densities, but constitute a 
very useful tool to investigate polarization among sea quarks. Specifically, 
we have found that SMC and Hermes data show a clear preference for par- 
ton distributions with 'up' antiquarks polarized along the nucleon spin and 
strange quarks with the opposite polarization Q . 

Regarding 'down' antiquarks, although the global fits attain marginally 
lower values of \ 2 when these are polarized in opposition to 'up' antiquarks, 
the differences are so small that we conclude that the data is unable to 
discriminate between the different alternatives. 

The agreement between polarized parton densities obtained from global 
fits and the constraints coming from hyperon decay data is particularly good 
in NLO, with departures of the order of a few percents, but larger for LO 
sets, partly due to higher order QCD effects. 
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Appendix: Parameters of the fits 



We present here the parameters of the best fits found in our analyses, 
both at LO and NLO. The parameters marked by an asterisk are the ones 
externally fixed from positivity constraints at large x. 



Parameter 


i + 


i - 


ii + 


ii - 


iii + 


iii - 




-0.021 


-0.021 


-0.034 


-0.030 


-0.014 


-0.012 


esu(3) 


-0.002 


-0.001 


0.056 


0.045 


-0.064 


-0.068 




0.873 


0.877 


0.921 


0.918 


0.936 


0.933 


0* 


3.20 


3.20 


3.20 


3.20 


3.20 


3.20 


In 


14.87 


14.66 


13.56 


13.69 


14.46 


14.40 


S u 


1.06 


1.06 


1.00 


1.00 


0.97 


0.97 


ad 


0.460 


0.460 


0.474 


0.465 


0.446 


0.443 


n 


4.05 


4.05 


4.05 


4.05 


4.05 


4.05 


id 


14.99 


14.99 


14.81 


14.99 


14.99 


14.99 


s d 


1.65 


1.67 


1.66 


1.66 


1.62 


1.63 




-0.062 


-0.062 


-0.069 


-0.069 


-0.063 


-0.063 




2.49 


2.49 


2.50 


2.50 


2.49 


2.49 


(is 


10.00 


10.00 


10.16 


10.16 


10.32 


10.32 




0.065 


0.046 


0.059 


0.044 


0.055 


0.042 


OLu 


1.01 


1.12 


1.25 


1.15 


1.12 


1.17 


N g 


0.25 


0.27 


0.40 


0.40 


0.7 


0.7 


a g 


1.5 


1.5 


1.5 


1.5 


1.5 


1.5 


ft, 


6.0 


6.0 


6.0 


6.0 


6.0 


6.0 


x 2 


330.52 


330.00 


330.98 


330.22 


332.31 


331.64 



Table Al: Coefficients for NLO total sets 
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Parameter 


i + 


i - 


ii + 


n - 


iii + 


m - 


esj 


-0.131 


-0.131 


-0.133 


-0.135 


-0.14 


-0.136 


e SU(3) 


0.071 


0.067 


0.032 


0.047 


0.036 


0.039 


Q-u 


0.768 


0.777 


0.766 


0.768 


0.77 


0.767 


P*u 


3.05 


3.05 


3.05 


3.05 


3.05 


3.05 


In 


14.88 


14.88 


14.97 


14.81 


14.98 


14.99 


Su 


1.64 


1.64 


1.69 


1.69 


1.69 


1.69 


ad 


0.5 


0.499 


0.502 


0.52 


0.532 


0.527 


ft 


3.77 


3.77 


3.77 


3.77 


3.77 


3.77 


Id 


14.00 


1 A OO 


10. u 


lo.z4 


lo. 1 


lo.yo 


5 d 


1.68 


1.75 


1.69 


1.69 


1.69 


1.69 


N s 


-0.069 


-0.07 


-0.068 


-0.070 


-0.069 


-0.069 


a s 


2.37 


2.4 


2.19 


2.25 


2.30 


2.32 


0s 


11.42 


11.42 


11.42 


11.42 


11.42 


11.47 




0.060 


0.043 


0.058 


0.042 


0.056 


0.041 


% 


1.19 


1.19 


1.2 


1.2 


1.2 


1.21 


N 9 


0.35 


0.35 


0.6 


0.6 


0.7 


0.7 


a g 


1.8 


1.78 


1.8 


1.8 


1.8 


1.8 


0* 


6.0 


6.0 


6.0 


6.0 


6.0 


6.0 


x 2 


336.24 


335.54 


334.19 


334.03 


334.26 


333.69 



Table A2: Coefficients for LO total sets 



26 





NLO 


LO 


Parameter 


1 


n 


iii 


i 


ii 


m 




-0.015 


-0.018 


-0.0025 


-0.126 


-0.126 


-0.126 


e SU(3) 


-0.003 


0.037 


-0.063 


0.016 


0.018 


0.009 


a u 


0.876 


0.915 


0.928 


0.766 


0.77 


0.771 


P*u 


3.2 


3.2 


3.2 


3.05 


3.05 


3.05 


In 


14.53 


13.37 


14.38 


14.99 


14.94 


14.84 


Su 


1.06 


1.0 


0.97 


1.66 


1.68 


1.68 


ad 


0.46 


0.47 


0.43 


0.471 


0.477 


0.49 


Pa 


4.05 


4.05 


4.05 


3.77 


3.77 


3.77 


id 


14.99 


14.85 


14.66 


15.0 


14.98 


14.97 


s d 


1.675 


1.66 


1.55 


1.72 


1.68 


1.68 


N s 


-0.062 


-0.069 


-0.064 


-0.065 


-0.068 


-0.069 


a s 


2.49 


2.5 


2.49 


2.38 


2.38 


2.38 


Ps 


10.0 


10.16 


10.32 


11.42 


11.42 


11.42 




0.25 


0.4 


0.7 


0.35 


0.59 


0.7 


OLg 


1.5 


1.5 


1.5 


1.8 


1.8 


1.75 


p; 


6.0 


6.0 


6.0 


6.0 


6.0 


6.0 


x 2 


228.66 


228.69 


230.33 


237.55 


236.1 


236.44 



Table A3: Coefficients for LO and NLO inclusive sets 
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